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ABSTRACT 

A new study of the diffuse Galactic 7-ray continuum radiation is presented, using a cosmic -ray propagation 
model which includes nucleons, antiprotons, electrons, positrons, and synchrotron radiation. Our treatment of the 
inverse Compton scattering includes the effect of anisotropic scattering in the Galactic interstellar radiation field 
(ISRF) and a new evaluation of the ISRF itself. Models based on locally measured electron and nucleon spectra 
and synchrotron constraints are consistent with 7-ray measurements in the 30-500 MeV range, but outside this 
range excesses are apparent. A harder nucleon spectrum is considered but fitting to 7-rays causes it to violate 
limits from positrons and antiprotons. A harder interstellar electron spectrum allows the 7-ray spectrum to be 
fitted above 1 GeV as well, and this can be further improved when combined with a modified nucleon spectrum 
which still respects the limits imposed by antiprotons and positrons. 

A large electron/inverse Compton halo is proposed which reproduces well the high-latitude variation of 7-ray 
emission; this is taken as support for the halo size for nucleons deduced from studies of cosmic-ray composition. 
Halo sizes in the range 4-10 kpc are favoured by both analyses. The halo contribution of Galactic emission to the 
high-latitude 7-ray intensity is large, with implications for the study of the diffuse extragalactic component and 
signatures of dark matter. The constraints provided by the radio synchrotron spectral index do not allow all of the 
<30 MeV 7-ray emission to be explained in terms of a steep electron spectrum unless this takes the form of a 
sharp upturn below 200 MeV. This leads us to prefer a source population as the origin of the excess low-energy 
7-rays, which can then be seen as a continuation of the hard X-ray continuum measured by OSSE, GINGA and 
RXTE. 

Subject headings: cosmic rays — diffusion — Galaxy: general — ISM: general — gamma rays: observations — 
gamma rays: theory 



1. INTRODUCTION 

Despite much effort the origin of the diffuse Galactic con- 
tinuum 7-ray emission is still subject to considerable uncer- 
tainties. While the main 7-ray production mechanisms are 
agreed to be inverse Compton (IC) scattering, 7r°-production, 
and bremsstrahlung, their individual contributions depend on 
many details such as interstellar electron and nucleon spec- 
tra, interstellar radiation and magnetic fields, gas distribution 
etc. At energies above ~ 1 GeV and below ~ 30 MeV the 
domi nant physical mech a nisms are yet to b e established (see, 
e.g., ftunteret al. 1997|, [Skibo et al. 1997|, P ohl & Esposito 



showed that the 7-ray excess could indeed be explained in terms 
of inverse Compton emission from a hard electron spectrum. 

The situation below several MeV is also unclear; Skibo et al. 
(1997 ) showed that the di ffuse flux measured by OSSE below 1 
MeV (Purcell et al. 1996) can be explained by bremsstrahlung 



only if there is a steep upturn in the electron spectrum at low en- 
ergies, but that this requires very large energy input into the in- 
terstellar medium. A model for the acceleration of low-energy 
electrons has been proposed by Schlickeiser (1997). An analy- 
sis of the emission in the 1-30 Me V range, based on the latest 
COMPTEL data, has been made by|MS99d, who found that the 
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The spectrum of Galactic 7-rays as measured by EGRET 
shows enhanced emission above 1 GeV in comparison with cal- 
culations based on lo cally measured proton and electron spectra 
( |Hunter et al.l997| ). Mori (1997) and Gralewicz et al. (1997) 
proposed a harder interstellar proton spectrum as a solution. 
This possibilit y has been teste d using cosmic -ray antiprotons 
and positrons ( VISR98 , MS99e ). Another explanation has been 



proposed by Porter & Protheroe (1997) and Pohl & Espos- 
ito (1998), who suggested that the average interstellar electron 
spectrum can be harder than that locally observed due to the 
spatially inhomogeneous source distribution and energy losses. 
Pohl & Esposito (1998) made detailed Monte Carlo simulations 
of the high-energy electron spectrum in the Galaxy taking into 
account the spatially inhomogeneous source distribution, and 



predicted intensities are significantly below the observations, 
and that a point-source component is probably necessary. Solv- 
ing these puzzles requires a systematic study including all rele- 
vant astrophysical data and a corresponding self-consistent ap- 
proach to be adopted. 

With this motivation a numerical method and correspond- 
ing computer code ('GALPROP') for the calculation of 
Galactic cosmic -ray propagation has be en developed (Strong 
& Moskalenko 1998, hereafter SM98). Primary and sec- 
ondary nucleons, primary and secondary electrons, secondary 
positrons and antiprotons, as well as 7-rays and synchrotron 
radiation are included. The basic spatial propagation mecha- 
nisms are diffusion and convection, while in momentum space 
energy loss and diffusive reacceleration are treated. Fragmenta- 
tion and energy losses are computed using realistic distributions 
for the interstellar gas and radiation fields. Our preliminary re- 
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suits w ere presented in Strong & Moskalenko (1997, hereafter 



SM97 ) and full results for protons, Helium, posit rons, and elec- 



trons in Moskalenko & Strong (1998a, hereafter |MS98a| ). The 
evaluation of the B/C and 10 Be/ 9 Be ratios, evaluation of diffu- 
sion/convection and reacceleration models , and full details of 
the numerical method are given in SM98 . Antiprotons have 
been evaluated in the co ntext of the 'hard interstellar nucleon 
spectrum' hypothesis in [VISR98 The effect of anisotropy on 
the inverse Compton scattering of cosmic -ray electrons in the 
Galactic radiatio n field is described in Moskalenko & Strong 
(2000, hereafter frflSO0| ). As an application of our model, the 
Green's functions for the propagation of positrons from dark- 
matter p article annihilations in the Galac tic halo have been eval- 



uated in Moskalenko & Strong (1999b) 



The r a tionale fo r our approach was given previously ( SM98 
MS984 |MSR98| , [Strong, Moskalenko, & Reimer 2000j 



the propagation region is bounded by R = Rh, z = ±Zh be- 
yond which free escape is assumed. We take Rh = 30 kpc. 
For a given Zh the diffusion coefficient as a function of mo- 
mentum and the reacceleration parameters are determined by 
B/C. Reacceleration provides a natural mechanism to reproduce 
the B/C ratio without an ad-hoc form for the diffusion coeffi- 
cient. The spatial diffusion coefficient is taken as (3Dq{p/ po) S ■ 
Our reacceleration treatment assumes a Kolmogorov spectrum 
with S = 1/3. For the case of reacceleration the momentum- 
space diffusion coefficient D pr is related to the spa tial coeffi- 
cient ( |Seo & Ptuskin 1994 perezinskii et al. 1990| ). The in- 
jection spectrum of nucleons is assumed to be a power law in 
momentum, dq(p)/dp cx p~ 7 for the injected particle density, 
if necessary with a break. 

The total magnetic field is assumed to have the form 



Bi. 



B Q e 



-(R-R e )/R B -\z\/z B 



(1) 



The values of the parameters (Bq, Rb, zb) are adjusted to 
match the 408 MHz synchrotron longitude and latitude dis- 
tributions. The interstellar hydrogen distribution uses HI and 
CO surveys and information on the ionized component; the 
Helium fraction of the gas is taken as 0.11 by number. En- 
ergy losses for electrons by ionization, Coulomb interactions, 
bremsstrahlung, inverse Compton, and synchrotron are in- 
cluded, and for nucleons by ionization and Coulomb interac- 
tions following Mannheim & Schlickeiser (1994). The distribu- 
tion of cosmic-ray sources is chosen to reproduce the cosmic- 
ray distribution determin ed by analysis of EGRET 7-ray data 
( [Strong & Mattox 1996 ). The source distribution adopted was 
described in SM98. It adequately reproduces the observed 7- 
ray based gradient, while being significantly flatter than the ob- 
served distribution of supernova remnants. 

The ISRF, which is used for calculation of the IC emission 
and electron energy losses, is based on stellar population mod- 
els and COBE results, plus the cosmic m icro wave background 
(CMB), more details are given in Sect. 2.1. IC scattering is 
treated usin g the fo rmalism for an anisotropic radiation field 
described in 



MS00. 



Briefly, the idea is to develop a model which simultaneously 
reproduces observational data of many kinds related to cosmic- 
ray origin and propagation: directly via measurements of nu- 
clei, electrons, and positrons, indirectly via 7-rays and syn- 
chrotron radiation. These data provide many independent con- 
straints on any model and our approach is able to take advantage 
of this since it aims to be consistent with many types of obser- 
vation. We emphasize also the use of realistic astrophysical 
input (e.g. for the gas distribution) as well as theoretical devel- 
opments (e.g. reacceleration). The code is sufficiently general 
that new physical effects can be introduced as required. We 
aim for a 'standard model' which can be improved with new 
astrophysical input and additional observational constraints. 

Comparing our approach with the model for EGRET data 
by Hunter et al. (1997), which used a spiral-arm model with 
cosmic -ray/gas coupling, we concentrate less on obtaining an 
exact fit to the angular distribution of 7-rays and more on the 
relation to cosmic -ray propagation theory and data. 

With this paper we complete the description of our model 
by describing the 7-ray calculation, and make a new derivation 
of the ISRF 1 . The 7-rays allow us to test some aspects of the 
model, such as halo size, whic h come from the previous work 
based on nucleon propagation ( SM98 ). We then use the com- 
plete model to try to answer the question: what changes to the 
'conventional' approach are required to fit the 7-ray data, and 
which are consistent with other constraints imposed by syn- 
chrotron, positrons, antiprotons, etc. ? Although no final an- 
swer is provided, we hope to have made a contribution to the 
solution. 

For interested users our model including software and result 
datasets is available in the public domain on the World Wide 
Web 2 . 

2. BASIC FEATURES OF THE GALPROP MODELS 

The GA LPROP models have been described in full detail 
elsewhere ( |SM98| ); here we just summarize briefly their basic 
features. 

The models are three dimensional with cylindrical symme- 
try in the Galaxy, and the basic coordinates are (R, z,p) where 
R is Galactocentric radius, z is the distance from the Galac- 
tic plane and p is the total particle momentum. In the models 

1 Since the IC scattering is one of the central points in our analysis, we feel that the derivation of the ISRF deserves a short description which we place in Sect. |2.l| , 
while more details will be given in a forthcoming paper. 

2 http://www.gamma. mpe-garching.mpg.de/~aws/aws. html 

3 While the propagation model uses cylindrically symmetric gas distributions this influences only the ionization energy losses, which affect protons below 1 GeV, 
so the lack of a full 3D treatment here has a negligible effect on the 7r° -decay emission. For the line-of-sight integral on the other hand, our use of the HI and H2 
annuli is important since it traces Galactic structure as seen from the solar position. 



Gas related 7-ray intensities are computed from the emissiv- 
ities as a function of (R, z, E y ) using the column densities of 
HI and H2 for Galactocentric an nuli based on 21 -cm and CO 
surveys ( Strong & Mattox 1996[ ) 3 . Our 7r°-decay calculation 



is given in pS98 al In addition our bremsstrahlung and syn- 
chrotron calculations are described in the present paper in Ap- 
pendices [}]; together with previous papers in this series this 
completes the full presentation of the details of our model. 

In our analysis we distinguish the following main cases: the 
'conventional' model which after propagation matches the ob- 
served electron and nucleon spectra, the 'hard nucleon spec- 
trum' model, and the 'hard electron spectrum' model. The 
'hard spectrum' models are chosen so that the calculated 7-ray 
spectrum matches the 7-ray EGRET data. 

2. 1 . Interstellar radiation field 

Since Mathis, Mezger, & Panagia (1983), Bloemen (1985), 
Cox & Mezger (1989), and Chi & Wolfendale (1991) no cal- 
culations of the large-scale Galactic ISRF have appeared in the 
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literature despite the considerable amount of new information 
now available especially from IRAS and COBE. These results 
reduce significantly the uncertainties in the calculation, espe- 
cially regarding the distribution of stars and the emission from 
dust. In view of the importance of the ISRF for 7-ray mod- 
els, a new calculation is justified. Moreover, we require the full 
ISRF as a function of (R, z,v), which is not available in the 
literature. Our ISRF calculation uses emissivities based on stel- 
lar populations and dust emission; as in the rest of the model, 
cylindrical symmetry is assumed. The dust and stellar compo- 
nents are stored separately in order to allow for their different 
source distributions in the anisotropic IC scattering calculation 
(MSOC ). Here we give only a brief summary of our ISRF calcu- 



lation; a fuller presentation will be given in a separate paper (in 
preparation). The resulting datasets are available at the address 
given in the Introduction. 

The infrared emissivities per atom of HI and H2 are based 
on COBE/DIRBE data from Sodrowski et al. (199 7), com bined 
with the distribution of HI and H 2 described in [SM98| . The 
spectral shape is based on the silicate, graphite and PAH syn- 
thetic spectrum using COBE data from Dwek et al. (1997). 

For the distribution of the old stellar disk component we use 
the model of Freudenreich (1998) based on the COBE/DIRBE 
few micron survey. This has an exponential disk with radial 
scale length of 2.6 kpc, a vertical cosh 2 (z) form with scale 
height of 0.346 kpc, and a central bar. We also use the Freuden- 
reich single-temperature (T = 3800 K) spectrum to compute 
the ISRF for 1-10 /im to calibrate the more extensive stellar 
population treatment. Since the Freudenreich model is based 
directly on COBE/DIRBE maps it should give an accurate ISRF 
at wavelengths of a few /im and serves as a reference datum for 
the more model-dependent shorter wavelength range. 

The stellar luminosity function is taken from Wainscoat et 
al. (1992). For each stellar class the local density and absolute 
magnitude in standard optical and near-infrared bands is given, 
and these are used to compute the local stellar emissivity by in- 
terpolation in wavelength. The z-scaleheight for each class and 
the spatial functions (disk, halo, rings, arms) given by Wain- 
scoat et al. (1992) then give the volume emissivity as a function 



of position and wavelength. All their main-sequence and AGB 
types were explicitly included. 

Absorption is based on the specific extinction per H atom 
given by Cardelli, Clayton, & Mathis (1989) and Mathis (1990). 
The albedo of dust particles is taken as 0.63 (M athis, Mezger, 
& Panagia 1983) and scattering is assumed to be sufficiently in 
the forward direction as not to affect the ISRF calculation too 



much. Again the gas model described in 3M98 is used 



The calculated R- and z-distributions of the total energy den- 
sity are shown in Fig. [j] in order to illustrate the ISRF distribu- 
tion in 3D. 

3. SUMMARY OF MODELS 

We consider 6 different models to illustrate the possible op- 
tions available. They differ mainly in their assumptions about 
the electron and nucleon spectra. The parameters of the mod- 
els and the main motivation for considering each one are sum- 
marized in Table 1. The electron and proton spectra and the 
synchrotron spectral index for all these models are shown in 
Figs. H| and Q. 

In model C ('conventional') the electron spectrum is ad- 
justed to agree with the locally measured one from 10 GeV to 
1 TeV and to satisfy the stringent synchrotron spectral index 
constraints. We show that the simple C model is inadequate for 
7-rays; the remaining models represent various possibilities for 
improvement. Model HN ('hard nucleon spectrum') uses the 
same electron spectrum as in model C, while the nucleon spec- 
trum is adjusted to fit the 7-ray emission above 1 GeV. This 
model is tested against antiproton and positron data. In model 
HE ('hard electron spectrum') the electron spectrum is adjusted 
to match the 7-ray emission above 1 GeV via IC emission, re- 
laxing the requirement of fitting the locally measured electrons 
above 10 GeV. Model HEMN has the same electron spectrum 
as the HE model but has a modified nucleon spectrum to ob- 
tain an improved fit to the 7-ray data. Model HELH ('large 
halo') is like the HEMN model but with 10 kpc halo height, 
to illustrate the possible influence on extragalactic background 
estimates. Finally, in model SE ('soft electron spectrum') a 
spectral upturn in the electron spectrum below 200 MeV is in- 
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Table 1 

Parameters and objectives of models. 











Injection index 






Model 8 - 


GALPROP 
code 


kpc 


Do 

2 — 1 

cm s 


electrons 


protons 


He 


Motivation/Comments 


C 


19-004508 


4 


6 X 10 28 


1.6/2.6 b 


2.25 


2.45 


matches local electron, nucleon data and syn- 
chrotron; consistent with p and e+ constraints 


HN f 


18-004432 


4 


3.5 x 10 28 


2.0/2.4 b 


1.7 


1.7 


matches high-energy 7-rays using hard nucleon 
spectrum; inconsistent with p and e+ constraints 


HE C 


19-004512 


4 


6 x 10 28 


1.7 


2.25 


2.45 


matches high-energy 7-rays using hard electron 
spectrum 


hemn 


19-004526 


4 


6 x 10 28 


1.8 


1.8/2.5 d 


1.8/2.5 d 


optimized to match high-energy 7-rays using 
hard electron spectrum and broken nucleon spec- 
trum; consistent with p and e+ constraints 


HELH 


19-010526 


10 


12 x 10 28 


1.8 


1.8/2.5 d 


1.8/2.5 d 


HEMN with large halo 


SE 


19-004606 


4 


6 x 10 28 


3.2/1.8 c 


2.25 


2.45 


matches low energy 7-rays using upturn in elec- 
tron spectrum 



Propagation parameters are given in 



SM98| (C, HE, HEMN models: 15-004500; HELH: 15-010500; HN: 15-004100). All models 
except SE and HN are with reacceleration (Alfven speed va = 20 km s _1 ). Do is the diffusion coefficient at 3 GV (5 GV for HN model). 
SE: S = 1/3, no reacceleration. 

^Electron injection index shown is below/above 10 GeV. 

c Nucleon spectrum normalization is 0.8 relative to model C. 

^Injection index shown is below/above 20 GeV/nucleon. 

°Electron injection index shown is below/above 200 MeV. 

^5 = —0.60/0.60 below/above 5 GV, no convection. 
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FIG. 2. — Electron spectra as obtained after propagation in our models compared with direct measurements. Data: Taira et al. (1993) (vertical lines), Golden et 
al. (1984, 1994) (shaded areas), Ferrando et al. (1996) (small diamonds), Barwick et al. (1998) (large diamonds). Left: Thin solid line: C model, dashes: HE model, 
thick solid line: HEMN model, dots: SE model. Right: Electron injection spectral indices 2.0-2.4, no reacceleration. 




FIG. 3. — Synchrotron spectral index for selected models. Measurements by different authors are shown by boxes. Data references: Webber, Simpson, & Cane 
(1980), Lawson et al. (1987), Roger et al. (1999), Broadbent, Haslam, & Osborne (1989), Platania et al. (1998) (gTey boxes and open box), Reich & Reich (1988), 
Davies, Watson, & Gutierrez (1996). Note that the error bar given by Webber, Simpson, & Cane (1980) is probably too small due to the difficulties of low-frequency 
radio measurements. Left: Solid line: model C, dashes: model HE, dash^dot: model HEMN, dash-3-dots: model SE. Right: Electron injection spectral indices 
2.0-2.4 (from bottom to top), no reacceleration; electron spectra as in Fig. (right). 



voked to reproduce the low-energy (<30 MeV) 7-ray emission 
without violating synchrotron constraints. 

Even given the particle injection spectra we still have the 
choice of halo size and whether to include reacceleration. We 
have used reacceleration models here except for the more ex- 
ploratory cases HN and SE. The propagation is obviously also 
subject to many uncertainties. The modelling of propagation 
can however simply be seen as a way to obtain a physically 
motivated set of particle spectra to be tested against 7-ray and 
other observations; in the end we test just the ambient electron 
and nucleon spectra against the data, independent of the phys- 
ical nature of their origin. In this sense our investigation does 



not depend on the details of the propagation models but still 
retains the constraints imposed by antiproton and positron data. 

4. SYNCHROTRON EMISSION 

Observations of synchrotron intensity and spectral index pro- 
vide essential and stringent constraints on the interstellar elec- 
tron spectrum and on our magnetic field model. For this reason 
we discuss it first, before considering the more complex subject 
of 7-rays. 

The synchrotron emission in 10 MHz - 10 GHz band con- 
strains the electron spectrum in th e ~ 1—10 GeV range (see e.g. 
Webber, Simpson, & Cane 1980). Out of the plane, free-free 
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FIG. 4. — Proton spectra as obtained after propagation in our models compared with IMAX data and published estimates of the interstellar spectrum. Solid line: 
using power-law injection spectrum (models C, HE), HashpH linp- with hrpatr in injection spectrum at 20 GeV (model HEMN), dotted line: hard nucleon spectrum 
(model HN^ Vertical hars- TMAX Hirprt measured values ( Mpnn pt al 1 Q9'/[l Evaluations nf thp intprstpllar sppptrnnv shadpH arpa- haspd r^n TMAY data H vTpnn et 
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al. 2000), connected filled squares: Webber & Potgieter (19H9) and Webber (1998), connected open diamonds: based on LEAF and 1MP-S ( ieo et al. 1991 ) 



absorption is only important below 10 MHz (e.g. Strong & 
Wolfendale 1978) and so can be neglected here. In particular 
the synchrotron spectral index (T oc v~P) provides information 
on the ambient electron spectral index 7 in this range (approx- 



imately given by (3 = 2 



but note that we perform the 



correct integration over our electron spectra after propagation). 

While there is considerable variation on the sky and scatter 
in the observations, and local variations due to loops and spurs, 
it is agreed that a general steepening with increasing frequency 
from (3 = 2.5 to (3 = 2.8 — 3 is present. Webber, Simpson, & 
Cane (1980) found f3 = 2.57 ± 0.03 for 10-100 MHz. Law- 
son et al. (1987) give values for 38^408 MHz between (3 = 2.5 
and 2.6 using drift-scan simulations wh ich lead to m ore reliable 
results than the original analyses (e.g. Sironi 1974: (3 ~ 2.4). 



A re cent reanalysis of a DRAO 22 MHz survey (R oger et al. 
1999J linds a rather uniform 22 - 408 MHz spectral index, with 
most of the emission falling in the range (3 — 2.40 — 2.55. 
Reich & Reich (1988) consider /3(408-1420 MHz)= 3.1 after 
taking into account thermal emission. Broadbent, Haslam, & 
Osborne (1989) find /3(408-5000 MHz) - 2.7 in the Galactic 
plane, using far IR data to model the thermal emission, but re- 
mark that 3.0 may be more appropriate for a full sky average 
(cf. the Reich & Reich value). Davies, Watson, & Gutierrez 
(1996) find an index range for 408-1420 MHz of f3 = 2.6-3.3 
for a high-latitude band, and state that 3.0 is a typical value. Re- 
cent n ew experiments give reliable spectral indices up to several 
GHz ( [Platania et al. 1998 ); they used a catalogue of HII regions 
to account for thermal emission. 

Fig. [| summarizes these estimates of the Galactic nonther- 
mal spectral index as a function of frequency. Since the elec- 
tron spectrum around 1 GeV is steepened both by energy losses 
and energy-dependent diffusion, we can conclude from the 
low-frequency (3 ~ 2.5 that the injection spectrum must have 
7 < 2.0. In fact our models require an injection 7 = 1.6 — 1.8 
to compensate the steepening and give reasonable agreement 
with the observed (3(v). Since all the models we will describe 
are chosen to have this injection index in the energy range pro- 



ducing radio synchrotron, they are all consistent with the syn- 
chrotron index constraints. 

The comparison with models also depends on the z- 
distribution of the magnetic field, since this affects how the 
spectral index is weighted with z, and will give larger indices 
for larger extents of B due the spectral steepening with z. Since 
the z-variation of B is unknown and otherwise plays a rather 
secondary role in our model we use our predicted (3 just for 
a representative intermediate Galactic direction (I = 60° , b = 
10°), which is taken as typical of the data with which we com- 
pare. The analysis is quite insensitive to the choice of direction. 

We evaluate synchrotron emission in more detail only for 
our model with hard electron and modified nucleon injection 
spectra (HEMN) since this is preferred from our 7-ray anal- 
ysis in the following sections. The values of the parameters 
adopted in equation (jl|), Bq = 6.I/1G, Rb = 10 kpc, zb = 2 
kpc, were found to reproduce sufficiently well the synchrotron 
index (Fig. [3]), and t he absolute magnitud e and profiles of the 
408 MHz emission ( Haslam et al. 1982 ) as shown in Fig. |[ 
The thermal c ontribution in the plane at this freque ncy is only 
about —15% ([Broadbent, Haslam, & Osborne 19891). A signif- 



icantly smaller field would give too low synchrotron intensites 
as well as a spectral index distribution which disagrees with the 
data, shifting the curve in the (3(v) plot to the left. Rb is con- 
strained by the longitude profile, and zb by the latitude profile 
of synchrotron emission. 

For comparison, Heiles (1996) gives Bq — 5/iG for the vol- 
ume and azimuthally averaged (uniform + random) field at 
the solar position based on pulsar rotation measures and syn- 
chrotron data. Vallee (1996) gives similar values. Our B value 
follows from the attempt to include 7-ray information on the 
electron spectrum throughout the Galaxy and is consistent with 
these other estimates. The radial distribution and magnitude of 
the magnetic field is also consistent with that used by Broadbent 
et al. (1990), as shown in Fig. ^. 

Our model cannot reproduce the asymmetries in latitude or 
fine details of the longitude distribution of synchrotron emis- 
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FIG. 6. — Magnetic field (Btot) distribution at z = as used in our models 
(solid line) compared to parametrization by Broadbent et al. (1990) (dashes). 



FIG. 7. — Gamma-ray energy spectrum of the inner Galaxy (300° < / < 
30°, | £> | < 5°) compared with calculations in C model (the injection spectrum 
chosen to fit local electron spectrum and synchrotron index). Curves show the 
contribution of TP hrernsstrahlnna and 7r Q -decay and the total. Data: EGRET 
Strnna MatmY 1 QQfth rnMP TRT (Strong et al. 1999), OSSE (I = 0, 25°: 
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sion and this is not our goal. An exact fit to the profiles, in- 
volving spiral structure as well as explicit modelling of random 
and non-random field components, as in Phillipps et al. (1981), 
Broadbent et al. (1990), Beuermann, Kanbach, & Berkhuijsen 
(1985), is not attempted here. 

5. GAMMA-RAY DATA 

For comparison of longitude and latitude profiles we use 
EGRET data from Cycle 1-4 in the form of standard counts 
and exposure maps in 10 energy ranges (bounded by 30, 50, 
70, 100, 150, 300, 500, 1000, 2000, 4000, and 10000 MeV). 
The systematic errors in the EGRET profiles due to calibration 



uncer tainty and time-depend ent sensitivity variations are up to 
13% (Sreekumar et al. 1998). These dominate over statistical 



errors in the EGRET profiles. 

The contribution from point sources was removed using the 
following procedure. Point-like 7-ray excesses were deter- 
mined in four energy regimes (30-100, 100-300, 300-1000, 
and >1000 MeV) using a likelihood method (M attox et al. 
1996). A detection threshold similar to EGRET source cata- 
logs was applied, and ~ 280 sources were selected by com- 
paring with sources from the 3rd EGRET catalog ( Hartman et 
al. 1999). Using either the catalog spectral indices, or a -2.0 
spectral index when no spectrum could be obtained, the flux of 
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each source was subsequently integrated for the standard en- 
ergy intervals. For the three brightest sources on the sky (Vela, 
Crab, Geminga pulsar) the flux determined in each energy in- 
terval was used directly. The simulated count distributions of 
the selected sources were subtracted from the summed count 
maps of the EGRET Cycle 1-4 data. 

The 7-ray sky maps computed in our models are convolved 
with the EGRET point-spread function generated for an E~ 2 
input spectrum (the convolution is insensitive to the exact form 
of the spectrum). For spectral comparison at low latitudes it 
is better to use results based on multicomponent fitting, which 
accounts for the angular resolution of the instrument; here we 
use the results of Strong & Mattox (1996), synthesizing the 
skymaps of Galactic emission from their model components 
and parameters. Cross-checks between this approach and the 
direct method show excellent agreement. At high latitudes, 
where the convolution has negligible effect, we generate spectra 
directly from the EGRET Cycle 1-4 data described above. 

For energies below 30 MeV only spectral data for the inner 
Galaxy (330° < I < 30°, |b| < 5°) are considered, COMP- 
TEL: Strong et al. (1999), OSSE: Kinzer, Purcell, & Kurfess 
(1999). The COMPTEL low-latitude spectrum is a recent im- 
proved analysis which is about a factor 2 above that given in 
Strong et al. (1997) and there remains some uncertainty as dis- 
cussed in Strong et al. (1999); however the difference has neg- 
ligible effect on our conclusions. The COMPTEL high-latitude 
spectra are from Bloemen et al. (1999), Kappadath (1998), and 
Weidenspointner et al. (1999). 

6. MODEL C (CONVENTIONAL MODEL) 

We start with a 'conventional' model which reproduces the 
local directly measured electron, proton, and Helium spectra 
above 10 GeV (where solar modulation is small) and which also 
satisfies the synchr otron co nstraints. The propagation parame- 
ters are taken from SM98 . This model has Zh = 4 kpc, reac- 
celeration with va = 20 km s _1 and a normalization chosen to 
best fit the local electron spectrum above 10 GeV. A break in the 
injection spectrum is required to fit both the synchrotron spec- 
trum and the directly measured electron spectrum; we adopted 
a steepening from —1.6 to —2.6 at 10 GeV. The local electron 
spectrum, the synchrotron spectral index, the local proton spec- 
trum, and the 7-ray spectrum of the inner Galaxy, are shown in 
Figs. [| H H and0, respectively. 

Model C is based entirely on non-7-ray data but still approx- 
imates the 7-ray data within a factor 3 over three decades of 
energy (10 MeV to 10 GeV). It also sati sfies the limits imposed 
by antiprotons a nd positr ons ( MSR98 , Moskalenko & Strong 
1998b, hereafter MS98b ), though our new calculation 4 shows 
some deficit of positrons below w 10 GeV (Figs. ^| and |l0[ see 
also discussion in Section 7.1). It still remains a useful first ap- 
proximation to serve as the basis for the developments which 
follow. 

The fit to the inner Galaxy 7-ray spectra is satisfactory from 
30 to 500 MeV but a large excess in the EGRET spectrum rel- 
ative to the predictions above 1 GeV is evident, as first pointed 
out by Hunter et al. (1997). Simple rescaling of either electron 
or nucleon spectra does not allow the agreement to be signfi- 
cantly improved. Harder nucleon or electron spectra are there- 
fore investigated below. 

The model also fails to account for the 7-ray intensities be- 



low 30 MeV as observed by COMPTEL and OSSE; attempting 
to account for this with a steeper electron spectrum immedi- 
ately violates the synchrotron constraints, unless the steepen- 
ing occurs at electron energies below a few hundred MeV, as 
discussed in Section 7.3 (model SE). To prove this important 
point, we consider a series of electron injection indices 2.0- 
2.4 in a model without reacceleration and diffusion coefficient 
index S = 1/3; this conveniently spans the range of reason- 
ably simple electron spectra. Fig. | (right) shows these elec- 
tron spectra, and Fig. |] (right) the corresponding synchrotron 
indices. The corresponding gamma-ray spectra are shown in 
Fig. ||. It is clear that an index 2.2-2.3 is required to fit the 7- 
rays, while this produces a synchrotron index about 0.8 which 
is substantially above that allowed by the data. Although this is 
for a particular propagation model, note that any combination 
of injection and propagation which fits the 7-rays will have the 
same problem, except for extreme models (e.g. our SE model). 

We emphasize that it is the synchrotron constraint on the 
electron index which forces us to this conclusion; in the ab- 
sence of this we would be free to adopt a uniformly steep elec- 
tron injection spectrum to obtain a fit to the low-energy 7-rays. 
This is the essential differ ence between present and earlier work 
(e.g. Strong 1996 , SM97) where the consequences of the syn- 
chrotron constraints were not fully appreciated. 

7. BEYOND THE CONVENTIONAL MODELS 
7.1. HN model (hard nucleon spectrum) 

One possibility to reproduce the 7-ray excess above 1 GeV 
is to invoke interstellar proton and Helium spectra which 
are harder than those directly o bserved in the heliosphere 
( Gralewicz et al. 1997 , Mori 1997 ). Spatial variations in the nu- 
cleon spectrum are quite possible over the Galaxy so that such 
an option is worth serious c onsideration. This model has been 
studied in detail in MSR98 in relation to antiprotrons, so that 
here we just summarize the results and also extend the calcu- 
lation to secondary positrons. The 7-ray spectrum of the in- 
ner Galaxy for a model with a hard nucleon injection spectrum 
chosen to match the 7-rays (no reacceleration, proton and He 
injection index = 1.7) is shown in Fig. |TTJ The corresponding 
propagated interstella r proton spectrum is shown in Fig. [I| 

As pointed out in MSR98 , the same nucleons which con- 
tribute to the GeV 7-ray emission through the decay of ir°- 
mesons produce also secondary antiprotons and positrons (on 
the same interstellar matter). The harder nucleon spectrum hy- 
pothesis, t herefore , can be t ested with measurements of CR p's 



and e+'s (VISR98, VIS98b). Above T p ~ few 10 GeV (for a 



power-law proton spectrum) the mean kinetic energy of parent 
protons is about 10 times larger than that of produced secondary 
p's, and roughly the same holds for 7-rays, so 10 GeV p's and 
7's both are produced by ~ 100 GeV nucleons. This relation is 
also valid for secondary positrons. Such tests are therefore well 
tuned, and sample the Galactic-scale properties of CR p and He 
rather than just the local region, independent of fluctuations due 
to local primary CR s ources. 

The conclusion of MSR98, based on the p/p ratio as mea- 
sured by Hof et al. (1996), was that antiprotons provide a sen- 
sitive test of the interstellar nucleon spectra, and that a hard 
nucleon spectrum overproduces p's at GeV energies. In Fig. ^ 
the predicted flux of antiprotons is compared with new absolute 
p fluxes above 3 GeV from the MASS91 experiment (Basini et 



4 Small differences from MS98b are due to use of our new ISRF and magnetic field model which modifies the energy losses, and use of = 4 kpc instead of 3 
kpc, which changes the propagation slightly. 



STRONG, MOSKALENKO, AND REIMER 






-2-1 1 2 3 A 5 

10 10 1 10 10 1 10 10 10 1 10 

Energy, MeV 



FIG. 8. — Gamma-ray energy spectrum of the inner Galaxy (300° < I < 30°, |6| < 5°) compared with calculations for electron injection spectral indices 
2.0-2.4 (from left to right, top to bottom) as shown in Fig. Bl Curves show the contribution of IC, bremsstrahlung, and 7r° -decay, and the total. Data: as in Fig. M. 
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FIG. 9. — Interstellar antiproton flux for models (MSR98) compared with 
local measurements. Solid line: C, dotted line: HN, dashed 
(direct measurements): triang les- MA SS91 experiment ( Basini et al. 199^| ), 
other points: see references inMSR98l 
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FIG. 10. — Interstellar positron spectra for models compared with data. 
Solid line: model C, dotted lin e- HN dashed line- HF.MN Data (direct mea- 
surements): HEAT experiment ([Barwick et al. 1998). 
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FIG . 11. — Gamma-ray data as in Fig. M compared with HN model (hard nucleon spectrum chosen to fit EGRET data above 1 GeV). 



al. 1999). While the agreement is good for the normal nucleon 
spectrum, the hard nucleon spectrum produces too many an- 
tiprotons by a factor of ~ 5, well outside the error bars of the 
three data points. We conclude that such a hard nucleon spec- 
trum is inconsistent with the antiproton data. 

Fig. [l(^ shows the interstellar positron spectrum for the con- 
ventional an d hard n ucleon spectra, where we used the formal- 
ism given in MS98a. The flux for the conventional case agrees 
with recent data ( Barwick et al. 1998 ) at high energies, where 
solar modulation is not so important. For the hard nucleon spec- 
trum the flux is higher than observed by factor ^4; this provides 
more evidence against a hard nucleon spectrum. However this 
test is less direct than p due to the difference in particle type, the 
large effect of energy losses, and the effect of solar modulation 
at lower positron energies. 



Taken together, the antiproton and positron data provide 
rather substantial evidence against the idea of explaining the 
> 1 GeV 7-ray excess with a hard nucleon spectrum. 

Note that in these tests we assume that only the protons have 
a local spectrum which is different from that on Galactic scales. 
We assume that the propagation parameters can still be derived 
from B/C, so that implicitly the heavier nuclei C, O are not af- 
fected. We could alternatively adopt a picture in which the C, 
O etc. also have a local spectrum different from Galactic scales, 
but then fitting the B/C ratio would imply S > 0.6 for the in- 
dex in the diffusion coefficient (nonreacceleration case) which 
is certainly problematic for high-energy anisotropy (see SM98 
and references therein) and larger than predicted by standard 
diffusion theory. Therefore we consider the only case worth 
testing at this stage is the one where only the protons are af- 
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FIG. 12. — Gamma-ray data as in Fig. [7] compared with HE model (electron 
injection index —1.7). 
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FIG. 13. — Gamma-ray data as in Fig. ^ compared with HEMN model 
(electron injection index —1.8, and modified nucleon spectrum). 



fected. In any case, pursuit of more complex options is beyond 
the scope of this paper. 

7.2. A harder interstellar electron spectrum (HE & HEMN 

models) 

An obvious way to improve the fit to the EGRET data above 
1 GeV is to adopt a harder interstellar electron injection spec- 
trum. Such models will not match the directly-observed elec- 
tron spectrum above 10 GeV, but this is not critical since the 
large energy losses in this region mean that large spatial fluc- 
tuations are expected ( Pohl & Esposito 1998 ). Hence we relax 
the constraint of consistency with the locally measured electron 
spectrum. 

For the HE model we adjust the electron injection index 
and the absolute electron flux to optimize the fit to the inner- 
Galaxy 7-ray spectrum; also the absolute nucleon intensity 
component) was reduced slightly (factor 0.8) within the limits 
allowed by the proton and Helium data. The inner-Galaxy 7- 
ray spectrum is shown in Fig. [j~2|. This model with its harder 
electron injection index (—1.7) satisfies the synchrotron con- 
straints and also leads to a better fit to the 7-rays above 1 GeV, 
but produces too few 7-rays at energies below 30 MeV by a 
factor 2-4. In this case the additional low energy 7-rays must 
be attributed to another component as discussed in Section (7j| 
in the context of SE model. 

Pohl & Esposito (1998) use an electron injection index 2.0 
with a Gaussian distribution of 0.2 for their 7-ray model. This 
has the effect of an upwards curvature which is equivalent to 
a harder spectrum similar to ours. Quantitatively, from their 
Fig. 4, the difference in effective index for their spiral arm 
model going from a single index to the index with dispersion 
gives in fact a difference 0.2 at high energies. So the dispersed 
index is equivalent to a single index, 2.0-0.2 = 1.8, which is 
similar to our 1.7 (HE model). 

Since the fit of HE model to the EGRET detailed spectral 
shape is still not very good above 1 GeV we can ask whether it 
can be improved by allowing more freedom in the nucleon spec- 
trum also (model HEMN). Some freedom is allowed since so- 
lar modulation affects direct measurements of nucleons below 



20 GeV, and the locally measured nucleon spectrum may not 
necessarily be representative of the average on Galactic scales 
either in spectrum or intensity due to details of Galactic struc- 
ture (e.g. spiral arms). Because of the hard electron spectrum 
the required modification to the nucleon spectrum is much less 
drastic than in model HN. By introducing an ad hoc flatten- 
ing of the nucleon spectrum below 20 GeV, a small steepen- 
ing above 20 GeV, and a suitable normalization, an improved 
match to the inner Galaxy EGRET spectrum is indeed possible 
(Fig. [l3]). The spectral parameters are given in Table 1. For the 
modified nucleon spectrum (Fig. [|) we must invoke departures 
from cylindrical symmetry so that the local value still agrees 
with direct measurements at the solar position. However this 
modification of the nucleon spectrum must be checked against 
the stringent constraints on the interstellar spectrum provided 
by antiprotons and positrons (as in model HN). The predictions 
of this model are shown in Figs. ^ and [l(J As expected the pre- 
dictions are larger than the conventional model but still within 
the antiproton and positron limits. 

So far this is the most promising model (at least for 7-rays 
>30 MeV) and hence we consider it further by testing the an- 
gular distribution of the emission. The synchrotron predictions 
for this model and comparison with data were presented in Sec- 
tion K Figs. [Ti] and [l| show the HEMN model latitude and 
longitude 7-ray distributions convolved with the EGRET point- 
spread function, compared to EGRET Phase 1-4 data. The sep- 
arate components are also shown. Since the isotropic compo- 
nent is here regarded as a free parameter it was adjusted in each 
energy range to give agreement with the high-latitude intensi- 
ties in the latitude plots, and the same value was used in the 
longitude plots. 

In latitude the agreement is always quite good, in longitude 
the maximum deviation is about 25% in the 100-150 MeV, 
150-300 MeV, and 4-10 MeV ranges, typically the agreement 
is better than 10%. We believe, this is satisfactory for a model 
which has not been optimized for the spatial fit in each energy 
range, but anyway these figures allow the reader to judge for 
himself. In any case, it should be remembered that unresolved 
point-source components and irregularities in the cosmic-ray 
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FIG. 14. — Latitude distribution of 7-rays (0° < I < 360°) for model HEMN (thick solid line). Separate components show the contribution of IC (dashes), 
bremsstrahlung (thin histogram), 7r°-decay (thick histogram), horizonal line: isotropic background. EGRET data: dotted line. 
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FIG. 17. — Gamma-ray data as in Fig. M compared with SE model (electn 
injection spectrum with upturn below 200~MeV). 



source distribution are expected to lead to deviations from our 
cylindrically symmetrical model at some level. 

In this model the contributions to the spectrum of the in- 
ner Galaxy from IC and 7r°-decay are about equal at 100 MeV 
and 6 GeV, 7r°-decay dominates between these energies, and 
bremsstrahlung produces ^10% of the total. The comparison 
shows that a model with large IC component can indeed repro- 
duce the data. A profile for 70-100 MeV enlarged to illustrate 
the high-latitude variation (Fig. |l6|) shows that this model also 
accounts very well for the observed emission; we regard this as 
support for the large IC halo concept. 

Turning to high energies, consider the latitude and longitude 
profiles for 4000-10000 MeV; the agreement shows that the 
adoption of a hard electron injection spectrum is a viable expla- 
nation for the >1 GeV excess. The latitude distribution here is 
not as wide as at low energies owing to the rapid energy losses 
of the electrons; both HN and HEMN models reproduce the 
observed spe ctrum, a nd latitude and longitude profiles almost 
equally well (MS98b), and hence it is difficult to discriminate 
between them on the basis of 7-rays alone. Independent tests 



however argue against HN as described in Section 7.1 



It is interesting to note that in fitting EGRET data Strong 
& Mattox (1996) found that the IC component had a harder 
spectrum than expected (see their Fig. 4), which was quite puz- 
zling at that time. Also the study of Chen, Dwyer, & Kaaret 
(1996) at high latitudes found a hard IC component. These re- 
sults can now be understood in the context of the HE or HEMN 
model; a renewed application of the fitting approach with the 
new models would be worthwhile and is intended for the fu- 
ture. All these results can be taken as adding support to the 
'hard electron spectrum' interpretation of the 7-ray results, and 
for the idea that the average interstellar electron spectrum is 
harder than that measured in the heliosphere. 

If this model is indeed correct it then implies that 
bremsstrahlung plays a rather minor role at all energies, con- 
trary to previous ideas, with IC and 7r°-decay accounting for 
^90% of the diffuse emission. 

Although we have introduced rather arbitrary modifications 
to both electron and nucleon spectra to better fit the 7-ray data, 



we note two recent indications that add support to our approach 
from independent studies. Baring et al. (1999) recently pre- 
sented models for shock acceleration in SNR which produce 
very flat electron spectra quite similar to what we require in 
the present case. Further study of the possible link between 
these spectra is in progress. A completely independent line 
of evidence for a low-energy flattening of the proton spec- 
trum has recently been presented by Lemoine, Vangioni-Flam, 
& Casse (1998), based on cosmic -ray produced light element 
abundances. If the proton and He spectra do differ from that lo- 
cally measured, this could of course also apply other primaries, 
but an investigation of this is beyond the scope of the present 
work, which focusses on 7-rays. 

7.3. A steeper spectrum of electrons at low energies (SE 
model) or a population of MeV point sources in the 
Galactic plane ? 

In order to reproduce the low-energy (< 30 MeV) 7-ray 
emission via diffuse processes it is necessary to invoke steepen- 
ing of the electron spectrum below about 200 MeV to compen- 
sate the increasing ionization losses. A steep slope continuing 
to higher energies would violate the synchrotron constraints on 
the spectral index, as discussed in Section Bl For illustration 
we show a non-reacceleration model with an injection index 
-3.2 below 200 MeV, and -1.8 above (SE model). This fits 
both the COMPTEL and EGRET data (Fig. |l7|) while remain- 
ing consistent with the synchrotron constraints (Fig. ||). The 
synchrotron index increase occurs at frequencies < 10 MHz, 
below the range where useful limits can be set (see S trong & 
Wolfendale 1978). Although the behaviour of the electron dif- 
fusion coefficient a t energies below 100 MeV is quite uncertain 
(Bieber et al. 1994) the propagated spectrum is here dominated 
by energy losses, which severely limit the electron range, so 
that this is not critical for our model. 

In this model 70% of the emission is bremsstrahlung and 
30% IC at 1 MeV. This is the only model in the present work 
which can reproduce the entire 7-ray spectrum. A possible 
mechanism for acceleration of low-energy electrons has been 
proposed by Schlickeiser (1997). However the adoption of 
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FIG. 18. — Energy spectrum of 7-rays from high Galactic latitudes (|6| > 70°, all longitudes) for HE 
EGRET total intensity from Cycle l^t data. COMPTEL data: high-latitude total intensity (open boxes: 



Weidenspointner et al. 19991. 



□d HELH model fright ShaHed areas: 
, diamonds: Kappadath 1998[ crosses: 



such a steep low-energy electron spectrum has problems associ- 
ated with the very large power input to the interstellar medium 
( Skibo et al. 1997 ), and is ad hoc with no independent support- 
ing evidence. Moreover the OSSE-GINGA 7-ray spectrum is 
steep er than E~ 2 below 500 keV (K inzer, Purcell, & Kurfess 
1999) which would require an even steeper electron injection 
spectrum than adopted here. It is more natural to consider that 
the COMPTEL excess is just a continuation of the same com- 
ponent producing the OSSE-GINGA spectrum. Most probably 
therefore the excess emission at low energies is produced by a 
population of sources such as supernova remnants, as has been 
proposed for the d iffuse hard X-ray emissio n from the plane ob- 
served by RXTE ( |Valinia & Marshall 1998t >, or X-ray transients 
i n their low state as s uggested for the OSSE diffuse hard X-rays 
(Lebrun et al. 1999). The contribution from point sources is 



then about 70% at 1 MeV, the rest being IC. Above 10 MeV 
the diffuse emission dominates. Yamasaki et al. (1997) esti- 
mate a 20% contribution from point sources to the hard X-ray 
plane emission, the rest being attributed to young electrons in 
SNR, but these are still localized 7-ray sources rather than truly 
diffuse emission. 

A model with a constant electron injection index -2.4 can 
also fit the low-energy 7-rays, but conflicts with the synchrotron 
index and fails to reproduce the high-ener gy 7-rays. W hile this 
has been a popular option in the past (e.g., Strong 1996 ), it can- 
not any longer be considered plausible. 

We note that another possible origin for the low-energy 7- 
rays, synchrotron radiation from ^100 TeV electrons, has been 
suggested by Porter & Protheroe (1997). 

8. HIGH LATITUDE 7-RAYS AND THE SIZE OF THE ELECTRON 

HALO 

Gamma rays provide a tracer of the electron halo via IC emis- 
sion. In considering the HEMN model we showed that the high- 
latitude variation of EGRET 7-rays is in good agreement with 
our large IC halo concept (Fig. Hq, Section 7.2). Indication for 



a large 7-ray halo was also found by Dixon et al. (1998) from 
analysis of EGRET data. Although we have used Zh — 4 kpc 



as the standard value in the present wo rk we now test larger 
values to derive limits. Studies of 10 Be ( SM98 ) gave the range 
Zh = 4 - 12 kpc for nucleons. Webber & Soutoul (1998) find 
z h = 2 - 4 kpc from 10 Be and 26 Al data. Ptuskin & Soutoul 
(1998) find Zh — 4.9^2 kpc. The high latitude 7-ray intensity 
increases with halo size due to IC emission (though much less 
than linearly due to electron energy losses), so that at least an 
upper limit on the halo size can be obtained. Fig. |l8| shows 
the 7-ray spectrum towards the Galactic poles for Zh = 4 kpc 
(HEMN model), and 10 kpc (HELH model). z h = 10 kpc 
is possible although the latitude profile for the 100-150 MeV 
range is then very broad and at the limit of consistency with 
EGRET data (Fig. |l9|). Further the isotropic component would 
have to approach zero above 300 MeV, so that this halo size can 
be considered an upper limit. 

If the halo size is 4-10 kpc as we argue, the contribution of 
Galactic emission to the total at high latitudes is larger than pre- 
viously considered likely and has consequences for the deriva- 
tion of the diffuse extragalactic emission (e.g., S reekumar et al. 
1998). An evaluation of the impact of our models on estimates 
of the extragalactic spectrum is beyond the scope of the present 
work. 

9. LUMINOSITY SPECTRUM OF OUR GALAXY 

So far some 90 extragalactic sources have been observed with 
the EGRET tel escope and several with COMPTEL (e.g., H art- 
man et al. 1997, 1999). Most of these sources are blazars. Such 
data usually serve as a basis for estimates of the extragalactic 
7-ray background radiation. However, the number of normal 
galaxies far exceeds that of active galaxies, it is therefore inter- 
esting to calculate the total diffuse continuum emission of our 
Galaxy as an example. Using cosmological evolution scenarios, 
this can then be used as the basis for estimates of the contribu- 
tion from normal galaxies to the extragalactic background. 

The luminosity spectrum of the diffuse emission from the 
Galaxy is shown in Fig. ^(j, based on models HEMN (4 kpc 
halo) and HELH (10 kpc halo). The total 7-ray diffuse lumi- 
nosity of the Galaxy above 1 MeV is Lq = 7.1 x 10 39 erg s _1 
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FIG. 20. — Gamma-ray luminosity spectrum of diffuse emission from the whole Galaxy using HEMN model (left), and HELH model (right). Total is shown as 
solid line. Separate components: IC (dashes), bremsstrahlung (dots), and 7r°-decay (dash-dot). 



for the HEMN model and 9.2 x 10 39 erg s" 1 for the HELH 
model. Above 100 MeV the values are 5.4 x 10 39 erg s _1 
and 6.3 x 10 39 erg s^ 1 , res pectively. These values are higher 
than previous estimates (e.g. |Bloemen, Blitz, & Hermsen 1984 : 
[1.6 - 3.2] x 10 39 erg s _1 above 100 MeV) due to our large 
halo and the fact that our model incorporates the EGRET GeV 
excess. 

10. CONCLUSIONS 

We have carried out a new study of the diffuse Galactic 7- 
ray continuum radiation using a cosmic-ray propagation model 
including nucleons, electrons, antiprotons, positrons, and syn- 
chrotron emission. 

We have shown that 'conventional' models based on locally 
measured cosmic -ray spectra are consistent with 7-ray mea- 



surements in the 30-500 MeV range, but outside this range 
excesses are apparent. A harder nucleon spectrum alone is 
considered but fitting to 7-rays causes it to violate limits from 
positrons and antiprotons. A harder interstellar electron spec- 
trum allows the 7-ray spectrum to be fitted also above 1 GeV, 
and this can be further improved when combined with a mod- 
ified nucleon spectrum which still respects the limits imposed 
by antiprotons and positrons. This is our preferred model, and 
it matches the EGRET 7-ray longitude and latitude profiles rea- 
sonably in each energy band. 

Such a model produces only 25-50% of the 1-30 MeV emis- 
sion by diffuse processes. The constraints provided by the syn- 
chrotron spectral index do not allow all of the <30 MeV 7-ray 
emission to be explained in terms of a steep electron spectrum 
unless this takes the form of a sharp upturn below 200 MeV. 
Therefore we prefer a source population as the origin of the ex- 
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cess low-energy 7-rays, which can then be seen as an extention 
of the hard X-ray continuum measured by OSSE, GINGA and 
RXTE. This is a quite natural scenario since it is very likely that 
the hard X-rays are indeed from unresolved sources, and the 
switchover from source-dominated to diffuse-dominated has to 
occur at some point; we propose here that it occurs at MeV 
energies. 

The large electron/IC halo suggested here reproduces well 
the high-latitude variation of 7-ray emission, which can be 
taken as support for the halo size for nucleons deduced from 



independent studies of cosmic -ray composition. Halo sizes in 
the range Zh = 4 — 10 kpc are favoured by both analyses. 

Our models suggest that bremsstrahlung plays a rather minor 
role, producing not more than ~10% of the Galactic emission 
at any energy. 

A part of this work was performed while IVM held a Na- 
tional Research Council/NASA GSFC Senior Research Asso- 
ciateship. 



APPENDIX 

A. SPECTRUM OF ELECTRON BREMSSTRAHLUNG IN THE ISM 

In order to calculate the electron bremsstrahlung spectrum in the interstellar medium, which includes neutral gas (hydrogen and 
Helium), hydrogen-like and Helium-like ions as well as the fully ionized medium, we use the works of Koch & Motz (1959), Gould 
(1969), and Blumenthal & Gould (1970). Our approach is similar to that used by Sacher & Schonfelder (1984), but differs in some 
details. Throughout this Section the units h = c = m e = 1 are used. 

The important parameter is the so-called screening factor defined as 

~ ' (Al) 



2707' 

where k is the energy and momentum of the emitted photon, 70 and 7 are the initial and final Lorentz factor of the electron in the 
collision. If 6 — ► 0, the distance of the high energy electron from the target atom is large compared to the atomic radius. In this 
case screening of the nucleus by the bound electrons is important. Otherwise, for the low-energy electron only the contribution of 
the nucleus is important, while at high energies the atomic electrons can be treated as unbound and can be taken into account as free 
charges. 

The cross section for electron-electron bremsstrahlung with one electron initially at rest approach es, at high energies (70, 7, k ^> 



1), the electron-proton bremsstrahlung cross section with the proton initially at rest (Gould 1969). Therefore the contribution of 
atomic electrons at high energies can be accounted for by a factor of (Z 2 + N) in place of Z £ in the formulas for the unshielded 
charge, where Z is the atomic number, and N is the number of the atomic electrons. In the present paper we treat free electrons 
in the ionized medium in the same way as protons, which is an approximation, but it provides reasonable accuracy for the range 
3-200 MeV where the bremsstrahlung contribution into the diffuse emission is most important. In any case the contribution from the 
ionized medium is of minor importance in comparison with that of the neutral gas. 

A.l. Low energies (0.01 < E kin < 0.07 MeV) 

This is the case of nonrelativistic nonscreened bremsstrahlung, A = S/(2afZ 1 ' 3 ) 1, In the Born approximation (2nZaf / flo, 
2nZaf I ' (3 <C 1) the production cross section is given by eq. 3BN(a) from Koch & Motz (1959) 

da _ l6Z 2 r 2 af ( p + p 



i^/^n^ln ^ , (A2) 
dk 3 fcpo \p -p J 

where a/ is the fine structure constant, po and p are initial and final momentum of the electron in the collision, (3q and j3 are initial 
and final velocity of the electron, and Je is the Elwert factor, 

_ - exp(-27rZa / /A))] 
JE /3[l-exp(-27rZa///3)] ' 



which is a correction for the cross section eq. QA2J) at nonrelativistic energies 



A.2. Intermediate energies (0.07 < -Ekin < 2 MeV) 



For the case of nonscreened bremsstrahlung (A ^> 1) the Born approximation cross section is given by ( Koch & Motz 1959 
eq. 3BN): 



= £f E z 2 r 2 a P ( 4 2707 — + P ° 1 £ Q7 | £ 7o ££ o 
dk ' E e / kp \ 3 ° pip 2 pi p 3 pp 

"8707 . ,2lh 2 +PoP 2 , k ( 7o7+Po 7o7 + P 2 . „, 7o7 



+L 



3 PoP P0P 6 ZpoP \ Po P P Po 



(A4) 



18 



DIFFUSE CONTINUUM GAMMA RAYS FROM THE GALAXY 



where 
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L = 2 In 
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1~P 
7o7 



PoP - 1 



k 



The factor £ is given by 



Z2 



3 (fe-B Wn )/96 



l-0.3e 



-k/c 



(A5) 



(A6) 



where b — 0.07 MeV, c = 0.33 MeV, and the expression in square brackets is a correction for the contribution of N atomic electrons, 
which is negligible at £^m ~ 0.1 MeV, but becomes as large as that of the protons at E^in ~ 2 MeV (see also general comments 
at the beginning of Appendix |a|). The second factor, in round brackets, is a correction to obtain a smooth connection between the 
approximations in the transition region near 0. 1 MeV and is essential only for small k. 

A.3. High energies (E^ n > 2 MeV) 
For the case of arbitrary screening we use eq. 3BS(b) from Koch & Motz (1959): 



da 
dk 



2 1 



1 + ^ 

7o 



If the scattering system is an unshielded charge, the functions are < 



6« = 4 



In 



-Hz, 

3 7o 
= Z 2 cj> l 
1 

" 2 



where 



(A7) 



(A8) 



For the case where the scattering system is a nucleus with bound electrons, the expressions for <j>i and <fi 2 are more complicated and 
depend on the atomic form factor. 

Corresponding expressions for one- and two-electron atoms (N = 1, 2) have been given by Gould (1969). Rearranging these one 
can obtain 



MN) 



(Z-N) 2 ^ 
(Z ~ N) 2 cj) u + 8Z 



8Zil- 



N-l 



Z 



1 R N (q) 2 
dq _ a (q - 8) 



1 - 



N — 1 



q° 



dq 



(A9) 



R N (q) 



where 



Ri(q) = l-Fi(q), 

R 2 (q)=2[l-F 2 (q)}-[1-F 2 (q)]/Z, 



Fl(q) = {l 
F 2 (q) = {1 



q 3 - 6S 2 q\n(q/S) + 36 2 q - AS 3 



q 2 /[2a f Z] 2 }- 2 - 
q 2 /[2a f (Z-5/16)} 2 }- 2 . 



(A10) 



Equations (|A9|), ( |A1C| ) are valid for any Z, including H~ ions. The formulas have been obtained under the assumption that the 
two-electron wave function of He-like atoms can be approximated by the product of one-electron functions in the form of Hylleraas 
or Hartree. For large A the expressions for <f>i and (f) 2 approach the unshielded value 

0i = 2 -> {Z 2 + N)(b u . 

For the case of neutral He atoms Gould (1969) gives also numerical values of <j>\ and <p 2 tabulated for the variable {6/ (2a/)} = 
0...10. The latter have been calculated for a Hartree-Fock wave function, which are considered to be more accurate than the Hylleraas 
function. At low energies (A ;> 2) both functions provide the identical results. 

At high energies where k, 70,7 3> 1, A ^ 4, eq. ( A4) with £ = 1 + N/Z 2 and = 1 can be applied, where electrons are treated 
as unbound in the same way as protons. 

A.4. Fano-Sauter Limit 

The formulas described above do not permit the evaluation of the cross section at the high-frequency limit k — * 70 — 1. The cross 
section obtained in the Born-approximation becomes ze ro in this limit, whi le the value is non-zero. The corresponding expression 
has been obtained by Fano in the Sauter approximation (Koch & Motz 1959): 



FS 



= ATrZ 3 a)r 2 e 
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1 
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A.5. Heavier atoms 

For the electron bremsstrahlung on neutral atoms heavier than He we use the Schiff formula ( fcoch & Motz 1959j , eq. 3BN(e)): 



d(T ~2 2 1 
dk =2Z W fc 



1 + X; - -— \ fin M (0) + 1 - - arctan b 
7o 3 7o/ V b 



(A12) 



.2. 

7o 



2 , ,2s 4(2 -b 2 ) , 8 2 



where 



b = 



z 1 / 3 



M(0) = 



1 



HU' '" v "' <5 2 (l + fc 2 )' 
B. SYNCHROTRON RADIATION 



For synchrotron emission we use the standard formula (see e.g. Ginzburg 1979). After averaging over the pitch angle for an 



isotropic ele ctron distribution, this gives the emissi vity e(z/, 7) of a single electron integrated over all directions relative to the field 
in the form ( 3hisellini, Guilbert, & Svensson 1988| ) 



0,7) = 4V3nr e m e cv B x 2 I K i/3 (x)K 1/3 (x) - -x [Kl /3 (x) - K 2 /3 (x)} 



(Bl) 



in units of (ergs s 1 Hz where v is the radiation frequency, 7 is the electron Lorentz factor, vb = eB/ (2ixm e c), B is the total 
magnetic field strength, x = vj (37 2 ^b), and K z is the modified Bessel function of order z. 
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